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I. INTRODUCTION 



The large number of strongly-coupled many-scales degrees of freedom in fully developed turbulent flows renders any 
analytical and computational approach impractical for very high Reynolds numbers. An alternative avenue is to look 
for a statistical description which by itself presents a difficult task. Only a few rigorous results, as for example the 
4/5-law have been derived from the underlying Navier-Stokes equation, employing some idealised assumptions 
like homogeneity, isotropy, etc. On the other hand, a multitude of empirical facts are known from the experimental 
analysis Q , giving rise to a huge playground for phenomenological modelling . 

One such class of data-driven models arc random multiplicative cascade processes (RMCP) j^, which are motivated 
by Richardsons energy cascade picture ^. They emphasise the observed multifractal nature of turbulent fields by 
assigning independent and identically distributed random weights qj to a nested hierarchy of incrtial range scales 
rj <lj = L/X' < L, confined by the integral scale L and the dissipation length scale rj. Amplitudes Ilj ~ qi ■ ■ ■ Qj of 
the energy-flux density field at intermediate scales Ij are then given as a multiplication of random weights, so that 
their respective moments 

(n-) = {qy = (^^y (1) 

indicate multiscaling. Note that interpretationally the energy flux field n(a;) is set equal to the energy-dissipation 
field e{x) and that experimentally Ilj is identified with the coarse grained surrogate energy-dissipation field = 
e{x)dx. As a consequence of the scaling two-point correlations 

also come with the same scaling Q] . It is not only the multifractal aspect of these models which is confirmed by the 
data [Q , more subtle observed quantities like multiplier distributions [§-0 and Markovian Kramers- Moyal coefficients 
PI have also been quantitatively reproduced by random multiplicative cascade processes |l5[ . But in spite of their 
success and inherent simplicity, random multiplicative cascade processes do have their drawbacks: at first there is no 
link to the fundamental Navier-Stokes equation and, secondly, no dynamics is included in these geometrical models 
as they are statically formulated in one-, sometimes in multi-dimensional spaces only. In this work we do not address 
the first problem, but present a solution for the second. 

A dynamical generalisation of the geometric random multiplicative cascade processes is proposed, which is contin- 
uous and causal in space and time, and which generates positive- valued multifractal fields, like the energy-dissipation 
field of fully developed turbulence. This paper is an extension of Ref. jl6| (see also [0). The model construction is 
presented in Sect. II and relies on the multifractal scaling form (|^) of two-point correlations only; stable statistics is 
used to facilitate analytical results, n-point correlations are discussed in Sect. Ill in the context of generalised fusion 
rules and are interpreted in terms of the random multiplicative weight picture of the geometric model precursors. 
Standard and non-standard obscrvables, like integral moments and multiplier distributions are addressed in Sect. IV 
and shown to agree with data. Some model variations are presented in Sect. V. 

II. MODEL-CONSTRUCTION BASED ON TWO-POINT CORRELATORS 

A. Ansatz 

For the positive- valued energy-dissipation field e{x,t), depending on one space and time coordinate x and t, we 
assume a spatio-temporal multiplicative process, 



f roc nOQ 

e{x,t)=expl dt' dx' f{x - x' ,t - t') -f{x' ,t') 

V-J — oo J —oo 



(3) 



where 7 ~ Sa{{dxdt)°' ^cr, — is distributed according to a stable white- noise field [|T8|. The normalisation 
{e{x,t)) = const ~ 1 determines the parameter fi — cr"/ cos(7ra/2). The symmetric index function 

1 iO< t<T, -g{T-t)<x<giT^t)) 
fix,t)^{ (4) 
(otherwise) 
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is defined by the kernel-function g{t); see Fig.l. Details of the explicit functional form for g{t) will be derived in Sect. 
II. B. The energy-dissipation field e{x, t) at spatio-temporal position (x, t) is composed multiplicatively of independent 
and identically distributed (i.i.d.) random-numbers exp {7(0;, t)}, lying inside an attached domain, bounded by the 
kernel-function g{t). The condition < t < T for the index- function / ^ ensures causality and defines the 
correlation-time T as the maximum of the temporal extension of the influence domain. L gives the maximum of 
the spatial extension and the correlation is given by 2g(0) = L. With the locality-condition g{T) = 0, the influence 
domain is attached to the point {x,t) in an unequivocal way. These boundary conditions for the kernel-function g{t) 
make the ansatz (|^) a dynamical, causal and continuous process for the energy-dissipation field. This process is also 
homogeneous in space and time, as we assumed an i.i.d random-field 7(2;, t) and independence of the form of the 
influence domain from the spatio-temporal position. 

The specification of the model is completed once the kernel- function g(t) and the remaining parameters of the 
random-field j{x, t) are determined. This will be done in the next subsection where multifractal statistics is demanded 
from the ansatz (y). Technically two-point correlations will be employed for this derivation. Their geometrical 
structure is illustrated in Fig. 2 for various spatio-temporal distances. We see that the correlation structure is a 
function of overlaps of infiuence domains, which on the other hand are bounded by the kernel-function g{t). Thus 
two-point correlations give direct information about the kernel-function. Since we want to construct our process to 
produce a multifractal field, we first have to specify the interplay between two-point correlations and multifractality. 

Multifractality of the energy-dissipation field e{x, t) in fully developed turbulence ||^ is usually observed by exam- 
ining the scaling behaviour {ei{x,t)") ^ Z""^" of the field amplitude ei{x,t) = l"^ J^^lj^ dx'e{x' ,t), coarse grained 
over a spatial domain of size / with centre x [Q; angular brackets denote the expectation value as the average over 
independent realizations of the field. Multifractality prevails once the scaling exponents t„ are represented by poly- 
nomials in n, although often for given n, it is already difficult to unambiguously extract a unique value for t„ via 
this standard route. Since the employed moments involve an integration over n-point statistics {e{xi,t) ■ ■ •e(a;„,t)), 
the latter are more fundamental to reveal multifractality. For example, an inertial range scaling (ef{x,t)'j ^ l^'^^ of 
the second-order integral moment automatically implies the same scaling {e{x^t)e{x + l,t)) ^ l^'^^ for the two-point 
correlation; however the reverse need not be true [ p^ , pO[ . 

For our purpose, we use two-point spatio-temporal correlators of arbitrary positive orders 

_ {e-^ix,t)e-^{x + Ax,t + At)) 
"^•"^^ ' ^" (e"i(x,t))(e"2(x + Aa;,i + At))- ^ ' 

Spatio-temporal multifractality is incorporated by demanding scaling relations for the temporal and spatial two-point 
correlators 

r„,,„,(Ax = 0,At) - At-^["i-"^l itL<At<T~tr,) (6) 

r„,,„,(Ax, At = 0) ~ Ax-^["i'"^l {v<Ax< xl) (7) 

with temporal and spatial multifractal scaling exponents T[ni,n2] and r [711,712]. For future purposes, we introduce 
arbitrary temporal At £ [tL,T — t,,] and spatial boundaries Ace £ [tItXl] for the scaling relations (^) and (0) to hold. 
One can think of [7;, xl] as defining the inertial range for two-point correlators. Since we are dealing with a dynamical 
model, [tL,T — t,,] then defines a temporal inertial range. 



B. Determination of the kernel-function g(t) 

The multiplicativity of the process with a stable white-noise field and the multifractal scaling of two-point 
correlators and are the only ingredients of our ansatz. The functional form of the kernel-function g{t) will 
now be deduced solely out of these premises. Fig. 2a illustrates the overlap of the influence domains of two points 
of the energy-dissipation field, having the same spatial position x but separated by a temporal distance At with 
II < At < T — Using the independence of 7(2;, t) at different spatio-temporal positions, the temporal two-point 
correlator (|^) becomes 

(exp I (Til + Jy(^^) dx'dt'-/{x', t') \ ) 

,XAx = o,At) = - ^— ^ ^ r . TY- 

exp |7ii Jyf^^t) dx'dt'-rix',t')jj (^exp jjia /y(At) dx'dt'j{x', t')jj 
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The influence domains of both points of the energy-dissipation field factorise into an overlap contribution V{At) and 
two non-overlap contributions V{0) — V{At). The latter appears in the denominator and the numerator of (^,) in the 
same respective orders ni and n2 and cancel out. The remaining contributions are those from the overlap of the two 
influence domains. The overlap V{At) is bounded by the kernel-function g{t). 
Using the Laplace-transform and the stable-property 

(exp {Ai7i + A272}) = cxp I ^^^^ (A? + A^) + /i (Ai + Asjj (9) 

with i.i.d. 7i,72 ~ Sa{a, ~1, ^) and Ai,A2 > 0, together with A7 Sa{\X\a; P,0) and 7 + /i ^ ^^(cr, /?, /i) for 
7 Saic, /3, 0), after a straightforward calculation we get for the temporal correlators 

r„„„,(Ax = 0, At) = exp i — ^ ((m + 7^2)" - < - <) ViAt) \ . (10) 

Eq. ( pO| ) is a purely geometrical relation connecting two-point statistics with overlapping domains. Since these overlap- 
domains are bounded by the kernel-function g{t), the scaling relation (^) translates into 

fT rait') 

V{At)=^ dt' dx'r^lnAt (11) 

JAt J~g(t') 



for ti, < At < T — tri- Combining ( |ll[ ) with (g) and (|lO[), the temporal scaling exponent then becomes T[ni,7i2] = 
a((ni -I- 712)" — 71" — 712), with some constant a. Finally after differentiating (pl| ) with respect to At, we get the 
kernel-function 

git) = ^ (12) 

Q, COS I } 

inside the temporal scaling regime t G [t^^T — t,^], where ao = — _2cr^ ■ The singularity of g{t) aX t ^ Q and the 
condition of locality g{T) = now justifles the introduction of the boundaries t^ and T — t,-i for the temporal scaling 
regime. 

The scaling relation (||) is valid independent of the behaviour of g{t) for t > T — tr,. This part of V[At) in Fig. 2a 
is always a complete part of the overlap as long as tL < At < T — t,; and therefore contributes only as a constant 
prefactor. This gives the possibility to achieve locality g{T) ~ with the arbitrary choice 

git) = - ^ (T - 1, <t<T), (13) 

t + t T + t ^ V 7 



where at t = T - t^, t = y'Tt^ - 7*^/4 - T + t^/2 is fixed by continuity with (|l|). For t real, we have to guarantee 

T > 7tri/4. The reason for the specific choice ( p^ is simply to stay in the same class of polynomials of order —1 and 
to be able to give an analytical expression for the new constant t. 

The temporal scaling behaviour is independent of the functional form of g{t) for < t < t^, as the intersection 
domain V{At) is not bounded by g(t) for t < t^. This domain is fixed by the equal time two-point correlators 

(cxp \{ni + 712) Jy,^. dx'dt'-/{x',t')\) 

r„„„,(Aa;,At = 0) = ^ ^ ' ' ^1 ^, (14) 

(exp |rii dx'dt'-f{x', t')jj {exp |n2 Jyi^^^) dx'dt'j{x',t')jj 

characterised by a spatial intersection domain V{Ax) and calculated in the same way as for temporal correlators; 
consult again Fig. 2b. For simplicity the boundaries rj = 2(7(r — tjj) and .t^ = '^gith) of the scaling regime (j^) are set 
equal to those of the temporal two-point correlators. Then the relation 

25(t(A.T)) = Ax (15) 

yields the intersection time t(A.T), where the boundaries of the two domains of influence intersect with each other. 
Making use of (112), the spatial overlap reads 
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V{Ax) = j dt'{2g{t') - Ax) 2 j dt'g{t') + 2ao f In - 1 j . (16) 



git) - 7^ (0 < t < tL) (17) 

t + tL 



The simple choice 



is not unique, but is acceptable as a simple continuous extension of (|1J) into the domain < t < t^, leading also 
to the perfect scaling relation (|^). Finally the cut-off is fixed by the condition g{0) = L/2. The kernel-function 
g{t) is completely specified by the scaling relations (|^) and (^. It comes with the four parameters T, L, II, tr}, with 
the first two representing the correlation time and length and the latter two limiting the perfect scaling range of the 
two-point correlators. Note that by construction, the limit of the temporal scaling regime 7^ is necessary for 
locality g{T) = and avoids ultraviolet divergence in the spatial correlations. Two more parameters are hidden in 
the field construction (^) and enter via the stable white-noise field; these are the stable index a and the parameter a, 
determining the absolute strength of the multifractal scaling exponents. By adjusting to the experimentally observed 
scaling exponents r[l, 1] = T2 = 0.225 Q, one of them, say a ~ 0.225/(2" — 2) becomes expressible in terms of the 
other. Fig. 1 shows the kernel-function g{t) for parameters a = 1.9 and T = L = 2^ , = trj = 10 in arbitrary units. 
Note again the three different functional forms for the kernel-function g{t), the arbitrary small-scale contribution for 
t > T ~ trj, the contribution II < t < T — t^j responsible for scaling of temporal two-point correlators and the large 
scale contributions t <tL responsible for the perfect scaling of spatial two-point correlators. 

Independent of a special choice of parameters, the presented solution for g{t) implies spatio-temporal equivalence 
of the scaling exponents 

T[ni, 712] = r[ni, 712] = a ((ni + 712)" - ~ rig ) . (18) 

This corresponds to the Taylor- Frozen-Flow Hypothesis 0| where correlations at a spatial distance are measured at a 
temporal distance. Here this hypothesis holds at least in the case of two-point correlations inside the scaling regime. 
Later on, it turns out that this is not true for arbitrary spatio-temporal n-point correlations. 



C. Link to RMCP 



In Ref. [|6| the kernel-function g{t) has been constructed in a different approach by linking the dynamical ansatz 
(^ directly to geometric random multiplicative cascade processes (RMCP). The functional form of the kernel-function 
g{t) in ||l^ is more or less the same, differing only in the regimes < t < and T ~ t,^ < t < T . 

In a discrete RMCP the energy-dissipation 

e{v) = \{q{h)=e^AY.\nq{l,)\ (19) 

resolved at the dissipation scale 77 is given as the product of i.i.d. multiplicative weights q{lj) of a nested hierarchy of 
scales 77 < Ij = L' /V < L' with A > 1. Comparing (|l|) with (|), the scale index j is replaced by a discrete time-index 
tj with 2g{tj) ~ Ij, leading to 

f rtj fg(t') 1 
q{lj)^cxp{ dt' dx'-f{x\t')\ . (20) 

[Jt,^, J-g{t') J 

Since 7(2;, t) is a stable white noise field, the independent multiplicative weights q(lj) become identically distributed 
once the volume associated to the time-strip integration, inside the exponential is scale-, i.e. j -independent. This 
leads directly to the result (|l|), vahd for = g~'^{L' /2) <t <T - t^. 
Identifying 



q{L) = exp |y dt' j 



9(t') 1 

dt'7(.T',i') (21) 

9(t') I 



and 
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q{ri) = expl dt'i dt'-f{x',t') 

[Jr-tr, J-g(t') 



(22) 



as large-scale and small-scale fluctuations, respectively, the field (||) can be written in the RMCP-analogue form 

where J = log;^ ((T — i^) /t^) represents the number of discrete cascade steps. 

Note that all q{lj) are i.i.d. for 1 < j < J, whereas q{ri) and q{L) come with different distributions. There 
is more analogy to geometric RMCP: in Ref. |2^, it has been demonstrated that the experimentally motivated 
unrestrictivc sampling of spatial two-point correlators leads to finite-size corrections to scaling, once the two-point 
distance approaches the integral length L and that this finite-size correction can be removed by a suitably tuned 
large-scale fluctuation. Since by construction, the dynamical RMCP (|^) comes with perfect scaling for the two-point 
correlators (|^), the large scale fluctuation q{L) has to be viewed as perfectly tuned. 

III. MORE ON N-POINT CORRELATIONS 

A. Spatio-temporal two-point correlations 

In the previous section we used only spatial and temporal correlations to fix the kernel-function g{t). Now we 
examine arbitrary spatio-temporal two-point correlations (e(x, t)e{x + Ax, t + At)). It is straightforward to calculate 
these correlations as a function of the spatio-temporal overlap V{Ax, At), which is shown as the shaded area in Fig. 
2c. Analytical results however need a proper distinction of several cases, which depend on the relative position of 
the two space-time points. All can be treated analytically. Instead of presenting all these formulas, we prefer only 
to illustrate them with Fig. 3. Fig. 3a shows the two-point correlations as a function of At and Ax. Parameters are 
a = 1.9, t[1, 1] = 0.225 and in arbitrary units tL = 5, trj = 1, T = L = 2^. For Ax = a rigorous straight- 
line behaviour is observed in this log- log-log plot for < At < T — t,,, which reflects the perfect scaling form 
of the limiting case (||). The same holds for the limiting case At = in the range rj < Ax < xl] see again 
Eq. (|). As Ax and/or At increase, more and more deviations from the rigorous scaling behaviour can be seen. 
Once Ax > L/2 + g{At) = Lcorr{At), the two influence domains do not overlap any more and as a consequence 
decorrelation {e{x,t)e{x + Ax > x + Lcorr{At),t + At)) = {e{x , t)) {e{x + Ax , t + At)) sets in. It is interesting to 
note that L^orriAt) displays the functional form of the kernel-function g{t). 

B. Equal-time n-point correlations I 

On our way to more complex correlation functions, equal time n-point correlations 

p„(xi,mi;...;x„,m„) = (e(xi,i)™i • • • e(x„, i)"'"} (24) 

of arbitrary order mi, . . . , to„ are now investigated. First, we will reveal their general structure. In the two following 
subsections they will be discussed in the context of fusion rules |^ and generic RMCP structures. We restrict 
ourselves to equal-time correlations. The case of temporal correlations is straightforward and does not give any new 
insight. The more complicated inspection of arbitrary spatio-temporal n-point correlation functions is briefly touched 
in Subsect. III.E. 

We again start with two points to exhibit the transition from two-point correlators to two-point correlation densities. 
Making use of (||) and the result (^, the equal-time two-point correlation density becomes: 



p(xi,mi;x2,m2) = (e(xi, <)"ie(.X2, t)™^) 

rm„m,{\x2 - x,\,At = 0) (e(xi,i)"^) (e(x2,0"=) 
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+ ((mi + 7712)" - mi - m2) (Vi n V2))} 



= p{m^\Vi/V2)p{m2\V2/Vi)p{m, + msl^i n V2). (25) 

Here K: = Jq dt 0^2^ represents the volume of the influence domain, attached to point i. Vi/Vj represents volume 

Vi without the overlapping volume ViCiVj resulting from Vj. In Fig. 4a the volumes V1/V2, V2/V1 and Vi fl V2 are 
labeled with their respective orders mi, m2 and mi + m2', these three disjunct contributions fill the union of the two 
influence domains with no remaining overlaps and consequently factorisc, leading to the last step of (p5|), where the 
abbreviation 

p{m\V) = exp I _ (m« - m) 1/| (26) 

has been introduced. 

For the third order n = 3 with points xi < X2 < we arrive at the straightforward generalisation of ( |25| ) 
p(a;i,mi;x2,m2;x3,m3) p(mi|Vi/V2) pim2\V2/{Vi U V3)) p{m3\V3/V2) 



p{mi + m2\{Vi n V2)/V3) -p[m2 + m3\{V2 n V7i)lVi) 

p(mi +m2 +m3|yi 0^2 nFg); (27) 

Consult Fig. 4b. For the most general case of an equal-time n-point correlation density with xi < X2 < . . . < a;„ we 
only state the result: 

n n—j+1 

p(xi,mi;...;a:„,m„) = Jl J| p(mfc + . . . + mfc+j_i|(Vfc n . . . n Vfc+j_i)/(Vfe_i U Vfc+j))- (28) 

j=i k=i 

One gets contributions of all common histories of neighbouring points, first the contributions every point has non- 
overlapping with the others, second all contributions from overlaps with one neighbour, with two neighbours,. . . and 
finally one contribution of the overlap, all points have in common. 



C. Equal-time n-point correlations II: generalised fusion rules 

After presenting the overall structure of equal-time n-point correlation densities, in terms of the intuitive picture 
with disjunct contributions from different overlap volumes of the respective influence cones, we now aim to rewrite 
them in terms of a fusion-rule picture. Again, the focus is flrst on two-point, then on three- and four-point, and finally 
on n-point correlation densities. 

Since the one-point moment (e(a;,i)™) is independent of the position x, the scaling behaviour of the equal-time 
two-point correlation density p{xi, mi; X2, m2) is identical to the two-point correlator (|^) for 77 < Ax = X2 — Xi < x^. 
Hence we arrive at 

p(a;i,mi;a;2,m2) ~ (a;2 - xi) ^'^I'^^l (29) 

with 

g[mi,m2] = 2a^!_2 ^^"^^ +"^2)° - - m^ ) . (30) 



For our purposes it is instructive to rcdcrive this result from the last line of (|25|). Noting that p{m\Vi/V2) 
p(m|l/i)/]o(m|X^i n V2) the scaling part of (E5|) is given by 



, . p(mi +m2|yi nV^) 

p{xi,mi;x2.,m2) —, . , , . - ■ 31) 

p(mi Vi n V2)p[m2\Vi n V2) 
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Realizing that Vi fl 1^ = V{/S.x = X2 — xi) ^ In Ax and making use of (p6|), Eq. ( ^T| ) then transforms into Eq. (29). 

In the expression ( p7| ) for the equal-time three-point correlator each of the six contributing factors is rewritten and 
then reordered, giving rise to 

p{xi,mi;x2,m2;x3,m3) = 

_ p{mi\Vi) 'p{m2\V2)p{m2\VinV3) pimslVs) 



pimi\Vi n V2) pim2\Vi n V2)pim2\V2 n V3) pim3\V2 n V3) 

p{rni + m2\Vi n V2) ^"^2 + m3\V2 n V3) _ ^ ^ 

= —, : 1,, ^-r,. =7 : „ -, . x P(mi + ma + m3\Vi H V3) 

p(mi -f TO2I Ki n K3) p[m2 + m3\Vi O I/3) 

p{mi + m2\Vi n V2) p(m-2 + "^31^2 n V3) p(mi -|- m2 -|- m3|Vi n V3)p(m2|Vi H ^3) 

p{mi\Vi n V2)p(to2|Vi n V2) ^(^21^2 n V3)p(to3|V2 n ^3) p(mi m2|l/i n F3)p(m2 + msl^L n T^s) 

^ {X2 - xi)-«t"^'™^l (X3 - X2)"«^'"^'"^l (X3 - . (32) 

The exponents £,[mi, 1712] and £,[m2, m3] are the same as in (^0[), whereas the third one is now given by i[mi, TO2, '713] = 
2^^ ((n^l + TO2 + ms)"" — ("^i + 1^2)°' — (to2 + TO3)" + wf" ). The last line of ( ^2| ) can be viewed as a generalised 
fusion rule; consult also Ref. 

The fusion rule (^) is easily generalised to arbitrary orders n. This is possible because all intersection domains 
can be written as a combination of ViHYj ^ In \ xj — Xi\, given that ry < \xj — Xi\ < xl for all pairs of points. Proven 
by complete induction the equal-time ri-point correlation density has the following structure: 

(n— 1 \ n— 1 n 

n - xo-^''"""'^^! n n - , (33) 
i=l J 3=2 l=j + l 

where 

C[mi, . . . , mi] = ^[toi + . . . + mi_i, mj] - ([m2 + . . . + mi_i,TOi] 

T2 



2" 



((toi -I- ... -I- iTii)" - (mi + . . . + mj_i)" - (to2 + . . . + ?7ii)" + (m2 + . . . + ?7i,_i)") (34) 



and Xl < X2 < . . . < x„ with r] < x^+i — Xi < xl- The equal-time n-point correlation function factorises into 
contributions at the smallest scales x^+i — Xi, next smaller scales Xi+2 — Xi ... and finally one contribution at the 
largest scale x„ — xi. The first product in ( p3| ) gives the contributions from the smallest scales involving the scaling- 
exponents ^ and the second product counts all other scales with modified scaling exponents ^. These modified 
scaling-exponents arise from the nested structure of the overlapping volumes. 



D. Equal-time n-point correlations III: RMCP interpretation 

So far the generalised fusion rule of the equal-time 71-point correlation densities has been expressed in two ways: 



Eq. ( |33| ) uses relative distances between the involved points, whereas Eq. (28) employs overlap volumes of respective 
influence domains. Since already in Section II. C we interpreted the influence domains in terms of random multiplicative 
weights, it is natural to interpret the generalised fusion rules also in terms of random multiplicative weights. In order 
to simplify the presentation, small- and large-scale fluctuations ( pT]) and (|2^) will be discarded. 

Using the assignment ([20|) the expression ([25|) for the two-point correlation density can be expressed as 



p(xi,mi;x2,m2) = (g™^'' {^'Y" . (35) 

In the RMCP language, 

na = (36) 

with 
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ftj rg{i') 

Vx= dt' / dx' = 2aolnA. (37) 
denoting the number of steps the two points have evolved together during the cascade history and 

. . ^1/(^10^2) V2/{Vinv,) 

n-.n- = (38) 

represent the number of cascade steps the two points go through independently after their branching. 
The three-point correlation density (E7^) displays the same structure: 



p{xi,mi;X2,m2]X3,m3) = {q""')^' (g™^)^'' (g™^')^'^ ^qmi+m^yi-^ ^^ms+ma^J^.a ^^m,+m2+msy 1,2.3 ^gg^ 

The first three factors constitute the multiplicative weights the three points do not have in common, the next two 
factors constitute the multiplicative weights which two of the three points have in common and finally the last line 
counts the multiplicative weights which all three points have in common. The numbers of multiplicative weights are 
given analogous to (^) and (|3^), 

. V1/V2 . V2/{Vi U V3) . V3/V2 

J, = = , J3 = (40) 

31,2 = y^ , J2,3 = (41) 

= y^ , (42) 

and arc related to each other by 

:h + jl.2 + jl,2,3 = 32 + .71,2 + j2,3 + .71,2,3 = h + .72,3 + jl,2,3, (43) 

reflecting the total number of RMCP cascade steps. Higher order correlation densities are straightforward and result 
in the same structure of common multiplicative weights of all combinations of neighbouring points. 



One more comment on the expression (39). In discrete RMCP, say with binary scale steps X ~ 2, three points 
first share a common cascade history with ji_2,3 steps. Then one point (say the third one) branches off and from 
then on evolves independently from the other two; consequently, j2,3 = could be zero and ja + ji,2,3 adds up to 
the total number of cascade steps. Only ji.2 could be different from zero, since the first two points still share some 
more common cascade history; once those two branches, ji — j2 independent steps are left until the RMCP evolution 
reaches the dissipation scale rj. However, this RMCP result, i.e. ji + ji,2 + Ji,2,3 = J2 + Ji,2 + ji,2,3 = js + ji,2,3, 
is not necessarily in conflict with the relation (|43|). The former reflects the ultrametric view of the nested hierarchy 
of RMCP length scales and is not observable for an experimentalist [T^- [r^ , |2l] , p3| , p^ . An unrestrictive sampling of 
discrete RMCP three-point correlations breaks the underlying ultrametric structure, leading to a modified result 
which may be closer to (^). This speculation should be tested with a discrete RMCP simulation. Definitely even 
more interesting would be to test the predicted generalised fusion rule structure (^), (|3^) and (^) of three-point 



correlations with data. Already the two-point correlations (0), (|25|), (U% and (35) are in excellent agreement with 



high Reynolds number data p5|, but three-point correlations could represent an even more stringent model test. 



E. Characteristic function 



In this subsection, we conclude our discussions on n-point statistics and examine the most general space-time 
correlation-function {e{xi,ti)^'^ ■ .. . ■ e{xn,tn)^"), containing all statistical information. It is equal to the generating 
function 

Z[^{x,t)] = (cxpl^J dx J dt ^{x,t) In e{x,t)^'^ (44) 

for the logarithmic field, once ^{x,t) = X]fe=i £,kS(t~tk)S{x — Xk) is chosen. The determination of Z is easily performed 
by regarding (H) as a summation of i.i.d. stable-distributions 7(0;, i) integrated over their influence domain: 
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Z[^{x,t)] = (cxpl^J dx' J dt'j{x',t')V^{x',t')^^ 



= exp —r— / dx' / dt' {V^{x', t') - V^{x\ t')) . (45) 

I V 2 j J —oo J — oc \ 

The first line is the resuh of inserting Eq. (^) into (^); also the abbreviation 

/>00 /"OO 

V^{x',t')^ dx dtf{x~x',t~t')^{x,t) (46) 



has been introdueed. For the second step of (^); the Laplace transform of stable distributions has been used. The 
result ( ^5|) is very convenient for numerical implementation: for ri-point correlation densities {e{xi, ti)^ . . . e{xn, ^n)^") 
the test function £,{x,t) — Y^^=i ^k5{x — Xk)5{t — tk) is a weighted sum of (5-functions, so that (Eq) simply counts 
the number of weighted spatio-temporal points inside the influence domain attached to (a;', tM. Note also, that the 
modified n-point correlation densities ((In e(a;i, ti))™i • • ■ (lne(a;„, i„) )"'"}, which follow from ( p4| ) by taking functional 
derivatives with respect to ^(x,t), only exist for the lowest orders, once a 7^ 2 is chosen for the random field 7(2;, t). 

IV. COARSE GRAINED OBSERVABLES 

In this Section, we examine statistical properties of coarse grained observables like integral moments and multiplier 
distributions. These are representative observables for the analysis of the experimentally measured surrogate energy 
dissipation field of fully developed turbulence. 

A. Integral moments 

The temporal integral moments of order n are defined as the expectation value of the coarse-grained temporal 
average of the energy-dissipation field 



ill i-to+t/2 

M(")(xo,to) = -My^ dt't{xo,t') 




(47) 



Within the inertial range of fully developed turbulence these moments show multifractal scaling with exponent Tt{n) 
Spatial integral moments, together with spatial scaling exponents Tx{n) are defined as the coarse grained spatial 
average of the energy-dissipation field. Invoking the Taylor-Frozen-Flow Hypothesis [Q, the spatial scaling exponents 
are believed to coincide with their temporal counterparts. For the time being, we present analytical model results for 
the temporal coarse graining and show numerical results for the spatial case. 

Integral moments of order n involve an integration over all n-point correlations. This integration can be done by 
integrating the second-order correlation function (^ with ni = 712 = 1: 

M?\xM=cj'-^\\c,m^ ' 



for < t < T — tff. The exponent T2 = a(2" — 2) should be set equal to the experimentally accepted value of 0.225 [Q. 
For tL <^ t <^ T ~ the last term ~ t^'^'^ dominates the second-order integral moment, so that {xq, to) ~ t^"^^ 
for intermediate scales. The analytic expressions for the coefficients c, ci and C2 are: 



2exp<{^y(T-i,)^ ( ^— ) (48) 



exp { 2^^(iL)} exp {^F(T - t,)} [T t,r 



cr = 2 ^ (2 - 2-^) ^-^ ) (49) 

exp{^y(t,.)} / ^2r. 4 X exp{^F(T-t„)}(T-t„)-tr . ^ ^ 

"^2 - 2 < — — — + 7; — — - 2 H - — > (50) 



l-2r2 V 2-2r2 2-2r2 / 2 - r2 
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with the V^(t)-expressions given by (^l]). The determination of higher order moments is straightforward, though 
cumbersome. 

Fig. 5 shows the local slope of logM^") as a function of the average domain logi and log/ for the temporal (Fig. 
5a) and the spatial integral moments (Fig. 5b), respectively. Parameters have been set a = 1.9 and in arbitrary units 
tL = 5trj = 5, T = L = 2*. For very small and very large scales there are weak deviations from perfect scaling, which 
increase as the order n increases. Comparing the temporal and spatial local slopes, there is clearly a larger scaling 
regime for the spatial moments. This is a result of the small scale statistics where temporal and spatial two-point 
correlations show a different behaviour. Once again, the message is that it is better to use two-point correlators 
instead of integral moments to extract multifractal scaling exponents. 



B. Multiplier-pdf 

Having examined moments of coarse-grained observables in the previous subsection, we now investigate the proba- 
bility density function p{M) of ratios of these observables, the so-called multipliers: 



Mt{x,tjt) 



(51) 

Clljle{x,t')dt'' 



Here we have only displayed the definition of a left-sided spatial and temporal multiplier Mx and Aft, respectively, 
where the daughter domain of size Z/2 is located at the left of the mother domain of size I. In the same manner right- 
shifted or centered multipliers, also with other scale steps, can be defined. A series of experimental investigations [p|jlC|] 
for large Reynolds number turbulent flows resulted in a /3-distribution with /3 ~ 3.2 for the left multiplier-distribution, 
scale-independent in the upper part of the inertial range. The most important result is reflected in the conditional 
multipHer distributions p{M {l)\M {21)) ||^, where correlations among multiphers with nested coarse graining domains 
at nearby scales have been observed: compared to the unconditioned pdf, for the left/right-sided multipliers one gets a 
more narrow pdf, if conditioned on a small parent multiplier and a more broader pdf, if conditioned on a large parent 
multiplier. On first view, the conditioning effects violate the assumptions of RMCPs since those require i.i.d. random- 
weights q{l). However, the contradiction is resolved once the observationally unavoidable unrestrictive sampling and 
the assumption of a nonconservative cascade generator is taken into account within the simple RMCPs . Since 

we pointed out that our dynamical approach can be viewed as a generalisation of these successful models, we now 
ask, for the outcome of these unconditioned and conditioned multiplier pdfs within the dynamical cascade model. 

We also report the differences between spatial and temporal multiplier distributions. The results presented here 
rely on simulations. Model parameters have been set r[l, 1] = 0.225, a = 1.92 and in arbitrary units T = L = 2^, 
tri = th = 1- The kernel-domain has been discretized into cells of size dt = 0.1, dx = 0.1, each filled with one 
random number 7. Stability of the results with respect to cell size and the number of independent samples of the 
energy-dissipation field has been checked. Fig. 6 shows (a) the temporal and (b) the spatial multipliers (|l]) obtained 
from 10^ field samples. For comparison the /3 distribution with (3 = 3.2 is depicted as the solid line. 

The unconditioned left-sided temporal multiplier-distribution is well-fitted by the experimentally observed (3- 
distribution, if the index of stability a is set equal to 1.92. For larger a all distributions become slightly more 
narrow and for smaller a they become slightly broader. Also the correct scale-correlations are reproduced in the 
conditioned distributions. These pdfs are scale invariant for It > that is over four binary orders of magnitude. 

For smaller scales all distributions become broader. The qualitative behaviour of temporal multiplier distributions 
remains unchanged once spatial multipliers are examined. But there is a clear quantitative difference, all spatial 
multiplier distributions are a little broader compared to their temporal counterpart. The differences arise again from 
the different small and large scale behaviour of n-point correlations. Scale independence and conditioning effects are 
nevertheless not affected. These model simulations demonstrate that the proposed spatio-temporal cascade process 
is also able to reproduce the experimental multiplier findings. 
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V. MODEL GENERALISATION TO 3 + 1 DIMENSIONS 



So far the dynamical cascade process has been formulated in 1 + 1 dimensions. Generalisations to n + 1 dimensions 
are straightforward. Their construction again relies on the scaling properties of spatial and temporal two-point 
correlators. In the following we only state the results for the special case n = 3 where the boundary conditions yield 
the most simplest analytical relations: 

Sticking to spherical symmetry, the kernel function g{t) of 1 + 1 dimensions is replaced by the time-dependent radius 
r{t), so that the scalar field f-{r,t) 



e(r, t) = exp ■ 



dt' 



It-T 



B'-", (r) 

T{t' -t + T)^ I 



dr'^{r' , t) 



(52) 



is then composed of time-dependent 3-dimensional balls -B^^]^ (r) . Scaling relations for spatial and temporal two-point 
correlators with symmetric scaling exponents T ~ t immediately imply 



r{t) 



^t+t T+t J 



1/3 
'iSL. 



T-t^<t<T 



tL<T-t^ 
Q<t<tL 



(53) 



The functional form for t > T — t^i again is arbitrary and only affects a constant prefactor in the temporal scaling. The 
functional form for < t <T — tr^ is fixed by the temporal scaling prerequisite y(Ai) In At. The new parameters 
^ and io are fixed by two boundary conditions, continuity at t = and a spatial scaling relation y(Ax) ~ In Ace: 



2S-2/3 



1-24 



3(1-20 



(54) 



The choice <^ = 2, for example results in ^ 0.33 and to ^ 0.5. Of course, other choices of ^ are also possible. 



VI. CONCLUSION 



The dynamical model presented here is constructed from a strict multifractal scaling of two-point correlators and 
employs stable statistics of the random-field 7. The latter allows an analytical treatment of n-point correlations. For 
fully developed turbulence, the dynamical model is not only able to reproduce the observed multifractal scaling of 
two-point correlators for the energy dissipation field, but also the observed multiplier distributions, including scale- 
correlations. The dynamic model also predicts specific generalised fusion rules for 71-point correlations; experimentally 
the sampling of at least three-point correlations should be feasible and could represent an even more stringent test of 
the proposed model. Another interesting observable, to be discussed in future, are 71-point cumulants of logarithmic 
field amplitudes [p4| . In this context, it might be necessary to formulate the proposed continuous dynamical model in 
a discrete way, introducing smallest spatio-temporal cells of size of the order of the dissipation scale and thus allowing 
to include deviations from log-stable statistics. 

Already in the last Section, we have sketched the formulation of the dynamical model in 3 -I- 1 dimensions. The 
extension to arbitrary n+1 dimensions is straightforward. Other generalisations in n+1 dimensions would be to include 
spatial anisotropy and boundary layer effects; work in this direction is already in progress p6|| . As possible applications 
of the proposed dynamical cascade model, we sec simulations for the spatio-temporal evolution of multifractal fields 
like rain and cloud fields in geophysics. In a numerically improved way, it might even serve as an efficient model 
generator for the turbulent small scale motion entering large-eddy simulations and other numerical approaches. 
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FIG. 3. Two-point correlations {e{x,t)e{x + Ax,t + At)) as a function of Aa; and At. Model parameters are quoted in the 
text. 
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FIG. 5. Local slopes of (a) the temporal and (b) the spatial integral moments A/'"' of order n = 2, 3, 4 as a function of 
the coarse-graining scale t and I, respectively. Model parameters are quoted in the text. For comparison the two-point scaling 
exponents r„ — 0.225(7i" — n)/{2" — 2) are indicated by straight lines. 
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